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Q\ ■ The Pseudo-Goldstone Boson (PGB) emission could provide a very efficient mechanism 
for explaining the cosmic Gamma Ray Bursts (GRBs). The PGBs could be produced 
during the merging of two compact objects like two neutron stars or neutron star - black 
hole, and then decay into electron-positron pairs and photons at distances of hundreds 
or thousands km, where the baryon density is low. In this way, a huge energy (up to 
more than 10 54 erg) can be deposited into the outer space in the form of ultrarelativistic 
e + e~ plasma, the so called fireball, which originates the observed gamma-ray bursts. The 
needed ranges for the PGB parameters are: mass of order MeV, coupling to nucleons 
9aN ~ few x 10~ 6 and to electrons g ae ~ few x 10 -9 . Interestingly enough, the range 
for coupling constants correspond to that of the invisible axion with the Peccei-Quinn 
symmetry breaking scale / ~ few x 10 5 GeV, but the mass of the PGB is many orders 
of magnitude larger than what such a scale would demand to an axion. Neither present 
experimental data nor astrophysical and cosmological arguments can exclude such an 
ultramassive axion, however the relevant parameters' window is within the reach of future 
experiments. Another exciting point is that our mechanism could explain the association 
of some GRBs with supernovae type Ib/c, as far as their progenitor stars have a radius 
~ 10 4 km. And finally, it also could help the supernova type II explosion: PGB emitted 
from the core of the collapsing star and decaying in the outer shells would deposit a 
kinetic energy of the order of 10 51 erg. In this way, emission of such an axion-like particle 
could provide an unique theoretical base for understanding the gamma ray bursts and 
supernova explosions. 
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1. Introduction 

The phenomenon of the so-called Gamma-Ray Bursts (GRBs) puzzles the theorists 
from many points of view (for a review see |TJ). The most striking feature of GRBs is 
that an enormous energy is released in few seconds in terms of gamma-rays having typical 
energies of few hundred keV. The energy emitted in GRBs is up to 10 53-54 erg, making 
it difficult to devise a mechanism able to transform efficiently enough the gravitational 
energy into such a powerful photon emission. In addition, the bursts show a variety of 
complex time-structures of the light curves. The time-structure of the prompt emission 
and the recent discovery of the afterglow fit the expectations of the so-called fireball model 
[0], in which the electromagnetic radiation is originated from the electron-positron plasma 
that expands at ultrarelativistic velocities undergoing internal and external shocks. The 
Lorentz factor of the plasma has to be very large, T ~ 100, in order to solve the so-called 
compactness problem || . So large values of T are difficult to achieve because they require 
a very efficient mechanism to accelerate the plasma. In addition, in order to avoid the 
contamination of the e + e~ plasma by more massive matter, the plasma has to be produced 
in a region of low baryonic density. 

Several mechanisms have been proposed to power the e + e~ plasma. In particular, 
it could be obtained via annihilation vv — > e + e~ of the neutrinos emitted by a heavy 
compact collapsing object - collapsar Alternatively, neutrinos can be produced at 

the merger of two compact objects, e.g. of two neutron stars or a neutron star and a 
black hole ||. The latter possibilities were recently analyzed in details in refs. P|-|rT||. Q 

The main challenge to the existing models is that the energy released in photons is 
astonishingly large. For example, the Gamma-Ray Burst GRB 990123 shows an energy 
release £ ~ 3.4 x 10 54 erg = 1.9M , if isotropic emission is assumed. Taking into account 
a possible beaming, this energy can be reduced down to £ ~ 6 x 10 52 erg There 
are also few other events (e.g. GRB 971214 and GRB 980703) with typical energies 
£ ~ 10 53 erg, which show no evidence for collimation ||14|| . The models invoking the 
vv — > e + e~ reaction as a source for the GRB have strong difficulties in reaching so large 
photon luminosities. Although during the collapse of compact objects a relevant amount 
of energy is normally emitted in terms of neutrinos, the low efficiency of the vv — > e + e~ 
reaction strongly reduces the energy deposited in the fireball. 

In the present paper we propose that light pseudoscalar particles - Pseudo-Goldstone 
Bosons (PGBs) - can be extremely efficient messengers for transferring the gravitational 
energy released in the merger of compact stars, into ultrarelativistic e + e~ plasma. We 
show that for certain parameter ranges the PGBs can be effectively produced inside the 
dense core of the collapsing system and then decay into 77 or e + e~ outside the system, 
in baryon free zones, thus giving rise to the ultrarelativistic e + e~ plasma. 

The advantage of assuming that the fireball is produced via the PGB decay instead of 
vv — > e + e~ annihilation is obvious. The latter process has a low efficiency - neutrinos 
deposit to plasma only few percent of the emitted energy and take the rest away. In 
addition, the process vv — > e + e~ can be effective only at small distances, less than 100 
km, which are still contaminated by baryon load, and it fails to provide a sufficiently large 



1 Recently some exotic mechanisms have been suggested related to the gravitational collapse of mirror 
stars 111. 
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Lorentz factor to the plasma, at most T ~ 5 |]9|,[T0[] . Instead, the PGB mechanism is 100 



percent efficient: decay can take place at distances about 1000 km and, in addition, all 
energy emitted in terms of PGBs is deposited to the e + e _ plasma with the Lorentz factor 
T ~ E/m e , where E is a typical energy of emitted PGBs (tens of MeV) and m e is the 



electron mass. This analysis is supported by the very recent result of ref. |15| , which 
shows that if the energy would be transfered to the plasma at distances of at least few 
hundreds km, a successful burst could be obtained with a Lorentz factor ~ 40 or so. 

The most familiar example of a PGB, the axion [16,17|, has already been considered in 
connection with GRBs |18| |. However, this was an invisible axion with the mass m a ~ 10~ 5 



eV. Such an axion has a lifetime much larger than the age of the universe and thus its 
decay is impossible. Instead, it was assumed that the axions could convert into photons in 
strong magnetic fields near the collapsing system. It is questionable, however, how large 
the magnetic field has to be in order to achieve a large efficiency for the conversion. 

In our analysis we phenomenologically consider a more general possibility, assuming 
that the PGB mass m a and its couplings to nucleons, electrons and photons are not 
constrained by the Peccei-Quinn relations. The paper is organized as follows. In section 
2 we shortly review the results of refs. PHn|] which describe the general scenario for 
mergers. In section 3 we consider the PGB emission as a source for the GRB and discuss 
the ranges for its mass and couplings which are compatible with the existing astrophysical 
constraints. In section 4 possible models for such axion-like particles are discussed. In 
section 5 we analyze possible signatures of the proposed mechanism. In particular, we 
show that the PGBs emission could generate a GRB associated with the supernovae type 
Ib,c, and also help the supernova type II explosion. Finally, in section 6 we summarize 
our findings. 



2. NS-NS and NS-BH mergers 



The scenarios depicted in refs. [0-11 1 for the NS-NS and NS-BH mergers are actually 
very similar. In both cases an axially symmetric structure develops. In the center the 
system collapses into a black hole of a few solar masses, while a fraction of matter (M ~ 
O.1M for NS-NS and M ~ O.5M for NS-BH) obtains enough angular momentum to 
resist immediate collapse into the black hole and remains in an accretion torus around the 
black-hole. Then this mass accretes from the torus into the BH going through a disk-like 
structure having a thickness which decreases approaching the BH. The central part of the 
disk is the most relevant for the radiation of energy, being the reservoir of largest density 
and temperature. The accretion rate from the torus dM/dt is of the order of 1 M Q /s for 
NS-NS and 5 M /s for NS-BH, so the duration of this phase is essentially similar in both 
cases, t ~ 0.1 s. 

According to refs. fl-ITH, the energy is radiated from the disk via the neutrino emis- 



sion. The deposition of energy in a region with low baryonic mass proceeds through vv 
annihilation into e + e~. The observed GRB is produced by the e + e~ plasma which ex- 
pands at ultra-relativistic velocity. The energy deposited into the e + e~ plasma ranges 
from E vV ~ 1Q48-49 er g j n ^ e case f NS-NS merger up to E uP ~ 10 516 erg for the NS-NH 
merger. The Lorentz factor are always rather small, T ~ 5. 

An important point in the simulations of refs. 0,0] concerns the density and temper- 
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ature profiles of the system. In ref. |J, where the neutrino trapping has been taken into 
account, the maximum density is ~ 10 11 g/cm 3 , with small regions reaching 10 12 g/cm 3 . 
On the contrary, in ref. |l(J where neutrinos are assumed to skip freely, densities are 
generally smaller by about an order of magnitude. When neutrino trapping is taken into 
account, larger densities, similar to the ones in ref. 0, are obtained. In ref. |L0| it is 
indeed shown that the density of the disk depends on the amount of energy lost by matter 
in this region. It is clear that only a totally self-consistent calculation can indicate if neu- 
trinos are trapped, reducing the luminosity of neutrino emission and therefore allowing 
for a larger density in the center of the disk. The connection between cooling and density 
is due to the contribution to the pressure coming from radiation and relativistic electrons 
and positrons. If the cooling is switched-off, only gas pressure is present and the system 
stabilizes at larger densities. 

It will be important in our discussion to keep in mind that a more efficient way of 
extracting energy from the system, e.g. via PGBs emission, will in turn modify the density 
in the center of the disk. Since in our work we do not pretend to solve the dynamics of 
the system, but we only indicate a new possibility for powering a strong GRB, we will 
consider for the central density of the disk values ranging from ~ 10 11 g/cm 3 down to 10 9 
g/cm 3 . 

The estimated temperature profiles of the disk are rather similar in refs. and [|K] 



reaching about 10 MeV in the central part and remaining above 1 MeV up to distances 
of the order of 100 km from the center of the system. 

3. GRB via PGB emission 

Let us consider the following Lagrangian describing the PGB couplings to the fermions 
ipi (i = e,p,n, ...) and electromagnetic field-strength tensor F^ v : 

£ = -Kn 2 a a 2 - ig ai aipa 5 ipi - -^-aF^F^ + .... (1) 

where in the last term the nucleon mass is taken as a regulator scale. In cases of 
the familiar models of the PGB like the axion, its mass and couplings are all determined, 
but for some coefficients, by only one free parameter /, the scale of the global U(1)pq 
symmetry. Instead, here we take a more phenomeno logical approach, considering m a , g ai , 
g a N and g ae as independent quantities. 

Let us start by discussing which constraints should be imposed immediately on the 
PGB mass and its couplings in order to obtain a strong GRB without violating the 
already existing experimental and astrophysical limits. In particular, we consider the 
mass and the coupling to the nucleon within the range 

0.3 MeV <m a < few MeV (2) 

10- 6 < g aN < 10~ 4 (3) 

The upper mass limit comes from the following. In order to have efficient PGB production 
in the dense systems under consideration, its mass should not over-exceed the character- 
istic temperature of the latter, typically of few MeV. The lower mass limit in (§) enables 
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to avoid the astrophysical constraints from the stellar evolution. The PGB production 
rate in the stellar cores with typical temperatures T up to 10 keV, is proportional to 
g 2 N exp(— m a /T). Thus, in order to avoid too fast stellar cooling, the exponential factor 
needs to be small if the constant g a ^ falls in the range indicated in ([$]). The lower limit 
on g Q 7v comes from supernovae. It implies that the PGBs are trapped in the collapsing 
core and thus the SN 1987A neutrino signal is not affected |nj.Q Finally, the upper limit 
in (Q) stands for a conservative interpretation of the set of constraints coming from the 
search of the decay K + — > n + a, from the beam dump and other terrestrial experiments 



20[ . The experimental limits on the couplings g ae and g ai will be discussed later. 

Let us discuss now in more details the range of the PGB parameters needed for suc- 
cessfully producing the GRBs. 

Mean lifetime 

The PGB can decay into photons and, if m a > 2m e , also into e + e~ . The corresponding 
decay width is r to t = r(a — > 77) + T(a — > e + e~), where 



A' 



a 2 ( 4m 2 \ 1/2 

r(a-e + e-) = ^m a . (5) 

8vr \ ml) 

As we anticipated in the introduction, we discuss the emission of PGBs by the central 
part of the toroidal system with typical size ~ 100 km, obtained during the merger of two 
compact objects. We assume that these PGBs take away a large fraction of the energy 
of the system and we want them to decay into photons or e + e~ preferentially outside the 
disk, in the region of low baryon density, thus giving rise to the hot relativistic plasma. 
Therefore a reasonable choice for the decay length of the PGBs, D a = ctj, is in the range 

100 km < D a < 10000 km. (6) 

Here r = is the lifetime of the PGB at rest and 7 = E/m a is its Lorentz factor. The 
average energy of the emitted PGBs can be taken as E ~ 3T, where T is a temperature 
of the system. Since T is of the order of several MeV, the Lorentz factor 7 ranges from 
2-3 up to 20 or more, depending on the mass of the PGB. The values of g ai and g ae 
corresponding to a given decay length D a , as functions of m a , are shown in Fig. 1. 

The upper limit in is not rigidly determined and is inferred from the following. In 
the literature it is argued that there is a correlation between some GRBs and supernovae 
type Ib/c |2l|j2^]. On the other hand, the observed frequency of supernova explosions is 



orders of magnitude larger than the observed frequency of GRBs. Therefore not every SN 
explosion is accompanied by a detected GRBs. If a supernova type Ib/c emits a GRB, 
the latter must be a rather weak one, so that it is detected only in few cases. Therefore, 
D a should be less than the typical radius of the type Ib/c supernova progenitors (R ~ 10 4 
km), so that the strength of the associated GRB is suppressed by a factor exp(—R/D a ). 



2 The region g a jy < 10 11 is also acceptable for the constraints coming from neutrino emission in super- 
novae jl9(| , but is not of interest for our analysis. 
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Figure 1. g ai and g ae as functions of m a [MeV], for the given decay length 
(D a = 10 2 , 10 3 and 10 4 km). The PGB energy taken is E = 15 MeV. 



Obviously, this suppression is enormous in the case of SN type II which progenitor has 
a typical radius larger than 10 7 km. In particular this makes clear the absence of a 
detectable GRB associated with SN 1987A. 



Mean free path 

The mean free path of a PGB having energy E depends on the density p and on the 
temperature T of matter, on the neutron and proton mass fractions X n and X p and, of 
course, on the PGB coupling to nucleons g a N- (Here and in the following, we take for 
simplicity the PGB couplings to proton and neutron to be equal, g ap = g an = g a N ). In 
ref. |23[ the mean free path is approximated by the following expression: 



\ a = g-*F(X n ,X p ,E/T) l a (p,T), 



(7) 



where the coefficient F 
energies E ~ 3T and 



1 + ZXnXp)- 1 (E/T) (1 + E/Ty 1 ' 2 is of order 1 for typical 



l a = (2.8 x 1(T 7 cm) x p- 2 T[MeV] 



1/2 



(8) 



with pi2 being the density in units of 10 12 g/cm 3 . 

For a given density and temperature one can estimate the critical value of g a N for 
which the PGBs become trapped. We denote this value as g f J N (p,T). Clearly, only an 
exact knowledge of the density and temperature profiles of the system can allow a precise 
determination of this quantity. However, we are only interested in a qualitative analysis 
and in our estimate we simply assume that the PGB is trapped when A a is smaller than 
the typical radius of the system R ~ 50 km. In other words, we define the critical value 
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g^ N as a function of a given density and temperature from the following equationP] 

(g t a r N r 2 l a (p,T) = R. (9) 

Hence, using eq. we obtain: 

g^ip, T) = 2.4 x lO-Vr^fMeV] 1 / 4 , (10) 

Assuming for instance T = 5MeV and p ~ 10 11 g/cm 3 , the PGBs are un-trapped for 
9aN < 9aN = 3-6 x 10~ 6 . If the density is smaller, g^ N increases. E.g., for p ~ 10 10 g/cm 3 , 
g ^ N = 3.6 x lfr 5 , while for p ~ 10 9 g/cm 3 , ^ = 3.6 x 10" 4 . 

Emission rate 

In the conditions of density and temperature typical for the central region of the disk, 
baryonic matter is non-degenerate, and the PGB production is dominated by bremsstrah- 
lung from nucleons. The analogous process from electrons is less efficient, also because 
the latter are partially degenerate at the densities and temperatures we are considering 
here. We use for the PGB production rate the expression [P~9|] : 

070^2 7-^/2 2 

Qnd = 9 2 aN R ND ( P ,T) = 9 2 aN ^% 9/ 2 = ^ 2 ^[MeV] 7 / 2 x 2.3-10 45 erg/cm 3 s(li; 

1UOTT / m ^ 

where, for simplicity, we consider symmetric nuclear matter, i.e. Kp n = Kp p = Kp, and 



a w ~ 15 is the pion-nucleon coupling constant JT9|. In Fig. 2 we show Log 10 (i?Ar£>[erg 
cm -3 s -1 ]) as a function of the temperature and of the density. 

Let us now estimate the maximal possible luminosity of the PGB emission as a function 
of g a N- Clearly, if the PGBs are emitted from a volume with given density and temperature 
in un-trapping regime, then their luminosity increases with g a ^ and reaches the maximum 
at g a N = Q^n^Pi T)- For the values of g a jy larger than g^ N , the PGBs become trapped and 
their luminosity starts to decrease, since only the PGB emission from the surface becomes 
possible. To compute the total luminosity we need also to estimate a corresponding volume 
associated with the un-trapped regime. 

The maximal luminosity for a volume V(p, T) with a given density and temperature 
reads therefore: 

LT X = i9aN(P, T)] 2 x R ND (p, T) x V(p, T) (12) 

We use the results of ref. to estimate the magnitude of the maximal luminosity. Let 
us consider, for example, the central zone with p ~ 10 11 g cm -3 and T ~ 3 MeV. The 
corresponding volume is roughly V ~ 10 6 km 3 . Then, according to eq. ( |T0D we obtain 
5a7v = 3-3 x 10~ 6 , and thus we get L™ ax ~ 10 56 erg/s, clearly a very large value which 
would be sufficient to power the most energetic bursts.^ It can be interesting to notice 

3 Taking into account that regions with high densities are smaller than R, the values we find for g*J N are 
clearly underestimated with respect to the actual values by a factor yf R/R(p, T) ~ few, where R(p, T) 
is the actual size of the region having a given density and temperature. However, this approximation is 
sufficient to our purposes. 

4 Certainly, such an enormous luminosity is unphysical. The actual luminosity should not exceed the rate 
of the mass transfer which, for NS-NS merger is dM/dt ~ 1 M© s -1 = 1.5 • 10 54 erg s _1 . As we noted 
above, in case of the NS-BH merger the later rate is about 5 times bigger. 
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Figure 2. Isocontours for Log 10 (-Rjvi)[erg/cm 3 s]) as a function of the temper- 
ature T[MeV] and of the density p[g/cm 3 ]. 



that Rnd is proportional to p 2 , while g f J N is roughly proportional to p _1 as is indicated 
by eqs. @ and (^). Hence, the maximal luminosity very weakly depends on the central 
density. 

Of course, the actual luminosity can be smaller. E.g., for g a ^ < g^ N the emission rate is 
simply proportional to g 2 N . If, on the other hand, g a N > <?aArG°> T), then the emission rate 
is suppressed due to the PGB trapping. In this case the volume in which the PGBs are 
trapped becomes ineffective, since now it contributes only through the surface emission, 
and so the dominant contribution comes from outer zones with smaller densities in which 
PGBs are still un-trapped. It is difficult to estimate the emission rate in the trapping 
regime. Taking, for example, g a N = 3 x 10 -5 , i.e. one order of magnitude larger then 
previous value, one can roughly estimate the total luminosity as L a ~ 10 55 erg/s, one 
order of magnitude less then above. 

Therefore, a luminosity larger than -^L™ ax ~ 5 x 10 54 erg/s can be obtained for g a ^ in 
the range: 

1(T 6 <g aN <3x 1(T 5 . (13) 

It is important to notice that, since the PGB luminosity is so large, a precise determination 
of the latter is inessential. Actually, for this range of g a N the PGB emission drains all the 
energy available from the system. 

For sake of simplicity in our analysis we have mainly considered the density and temper- 
ature profiles which are characteristic for the NS-NS merger. Let us shortly comment also 
the case of NS-BH merger. The latter gives somewhat larger densities and temperatures 
than the NS-NS merger. As observed above, a change in the density is not affecting the 
maximal luminosity L™ ax . A slightly larger temperature, on the other hand, corresponds 
to a larger production rate Qnd and therefore to a larger luminosity. 
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4. The PGB models 

Let us discuss now which particle physics candidate could be the PGB with the above 
properties. As we told already, the most familiar example is the axion, the PGB associated 
with the spontaneous breaking of the Peccei-Quinn symmetry U(1)pq needed for the 



solution of the strong CP-problem. The standard axion [16| is excluded by laboratory 



experiments long ago. However, several types of invisible axions have been considered 



in the literature: the DFSZ axion ||24j| , the KSVZ or hadronic axion ||25|| , archion p6 
etc. In these models the PQ symmetry is spontaneously broken due to the large VEV 
of some scalar S, singlet of the standard model: (S) = f/y/2. The axion interactions to 
matter fields have the form ([!]) where the coupling constants to the fermions and photon 
are related to the scale / as follows: 

mi a m N 

g a i = Ci— (i = e,p,n...), 9ay = c^— — , (14) 

with q and c 7 being the model dependent coefficients (generically, axion models employ 
a number of heavy fermion states that contribute to the colour N and electromagnetic 
E anomalies of the U(1)pq current). For any type of axion, its couplings to nucleons do 
not vary strongly from model to model and, typically, c P) „ are in the range 0.3-1.5. As 
for g ae and g a7 , their model dependence is stronger. For example, for the DFSZ axion 
model containing only the standard fermion families (N = Nf is the number of families, 
and E/N = 8/3) we have 

cf FSZ = si^p, c dfsz = ^ (15) 

1 ~~T~ Z 

where tan/3 = v u /vd is the VEV ratio of the up and down Higgs doublets H u and Ha- 
Thus, for Nf = 3 and a natural range v u /vd > 1, c e varies from 1/2 to 1, while c 7 
depends on the value of up/down quark mass ratio z = m u /md = 0.3 — 0.7. Therefore, 
Gael ' 9a~j = 0.47c e /c 7 can vary between 0.10 and 0.35. On the other hand, for the DFSZ-like 
axion we have g aN = (c N m N /c e m e )g ae ~ (10 3 - 10 4 )5f ae . 

For the hadronic axion or the archion the coupling to electrons is strongly suppressed. 
For example, for the hadronic axion g ea = at tree level and it emerges from radiative 



corrections. Namely, we have |I7| 




c y = E- 2{A + Z ] n, (16) 
7 3(1 + z) ' K ' 

where A ~ 1 GeV is a QCD scale. In the GUT context one typically has E/N = 
8/3, however this value is not mandatory for a general case. Essentially the same order 
estimates hold for the archion couplings, which in addition can have also a small tree level 
coupling to electron, with c e = 0{m e /m T ) ~ 10~ 3 [ ^B] . 



On the other hand, the axion mass is also related to the PQ symmetry breaking scale 
as m a ~ Nm-nfn/f and thus for / ~ 10 5 GeV it is too light to be applicable for our 
mechanism. However, we do not constrain the magnitude of m a with this relation and 
leave it as free parameter, in the range indicated in (fj). In the context of a theory 
having the U(1)pq symmetry, this would mean that the axion mass is not determined by 
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Figure 3. Possible values of g ae and g ai for a given decay length and axion 
energy: D a = 10 4 km, £=15 MeV (a) and D a = 500 km, E=9 MeV (b), for 
different values of the mass m a =0.5,l, 4 MeV. Shaded areas indicate typical 
correlations between g ae and g ai for the DFSZ axion (dark) and hadronic axion 
or archion (light) models. 



the colour anomaly but rather by some other dynamics. For example, one can imagine 
that its mass emerges from the U(1)pq current anomaly related to some hidden gauge 
sector with the confinement scale larger than A (in this case this axion would solve the 
strong CP problem in this sector). Alternatively, the Lagrangian could contain the higher 
order terms cutoff by the Planck scale which explicitly break the PQ symmetry and thus 



produce the PGB mass []27[]. In either case, such PGB would not be anymore relevant for 
the strong CP problem. One can rather consider it as a "would be" axion which could 
work for strong CP but is prevented to do so by the explicit U(1)pq breaking terms which 
in general have no reason to respect the strong 9 phase orientation. 

For example, taking the Planck scale induced terms in the Lagrangian as X(S + S) 2 S/Mp, 



where A is a coupling constant, one can estimate the PGB mass as [27 



m » = ^(lo4v) 3/!x6MeV (17) 

Although the coincidence most probably is simply accidental, we remark that for the 
DFSZ-like "failed" model with / = few x 10 5 GeV (and hence g ae ~ few x 10~ 9 , g a N ~ 
few x 10~ 6 ), the PGB mass in eq. flTT| ) naturally falls in the range of few Me Vs. 

We come now back to the analysis of our constraints on the PGB couplings, interpreting 
these constraints in the light of the above discussed models. In Fig. 3 we show the 
isocurves corresponding to a fixed decay length D a of the PGB with a given energy E, 
for different values of m a . We also indicate typical correlations between the couplings g ae 
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and g ai as they emerge for various axion models. 

We see that there can be two type of solutions. For the DFSZ axion it is definitely 
the case of decay into e + e~ which can work. It corresponds to values of the constant 
g ae ~ fewlO -9 , i.e. / ~ few x 10 5 GeV and axion mass of few MeV. Interestingly, the 



parameters range excluded by the present reactor experiments |28[ is approaching the 
parameters spot we are indicating, but the parameters value we are suggesting is not yet 
excluded. It must also be remarked that for such values of / the axion-nucleon coupling 
constant lies within the range of interest for the GRB, g a ^ ~ few x 10 -6 . 

Another solution for the DFSZ, via axion decay into photons (vertical lines in Fig. 3 
for m a < 1 MeV) leads to g ae > 10~ 7 which is ruled out by experiment 

As for the case of the hadronic axion or archion, which have a coupling to electrons 
suppressed by 3-4 orders of magnitude if compared to the DFSZ axion, the decay rate into 
photons is comparable to that into e + e~ when the latter is allowed (m a > 2m e ). However, 



it requires a too large g ai , above 10 7 . The experimental limit of ref. p9| |, / > 1 TeV 



excludes this possibility. Perhaps, this solution is still allowed with a larger axion mass, 
m a > 5 MeV, but this question needs an additional investigation. 



5. GRB and Supernovae 

In the previous sections we have considered GRBs emerging by axion emission during a 
NS-NS or NS-BH merger. However, the proposed mechanism can have strong implications 
also for the GRBs associated with supernova explosion, as well as for the supernova 
explosion itself. 

In the typical situation we considered, axions produced in a compact disk fly away 
taking the big portion of the available energy. In case of merger the typical size of the 
system is of the order of 100 km - there is essentially no baryonic matter outside this 
radius, so the axion flux, after decay into photons or e + e~ at a distance larger then 100 
km, converts into a fireball - a bubble (or jet) of hot relativistic plasma expanding with 
a large Lorentz factor.^ 

However, the axions can be produced also in the core of collapsing stars. In particular, 
for g a N ~ 10 -6 the axions are essentially in the trapping regime in the supernovae ex- 
ploding via core collapse. Axions are emitted from the axiosphere, having a temperature 
decreasing with time from T ~ 4 MeV at t = 100 ms to T ~ 3 MeV at t — 1 s. The 
corresponding luminosities are L a ~ 3 ■ 10 51 erg/s at t ~ 100 ms and L a ~ 0.5 ■ 10 51 
erg/s at t ~ 1 s ||23|| . Therefore, in total an energy ~ 10 51 erg can be emitted during 
the collapse period and subsequent cooling of the proto-neutron star, in terms of axions 
having an average energy E ~ 10 MeV. Since E is of the same order as the one in case 
of mergers, the value D a = cr r y is again of the order of few hundred or few thousand km. 
Then in our scenario the impact of the axions simply depends on the geometrical size of 
the collapsing star. 

5 A possible speculation concerns the observed spectrum of GRBs. The latter display a characteristic 
paucity for energies below few tens of keV. Moreover it typically has a maximum (the hardness H |30[ ) 
for energies of the order of few hundreds keV. Although it is not clear if there is any relation between the 
"internal engine" and the observed spectrum, we cannot avoid to remark a similarity between the typical 
value of the hardness and the mass of the PGB we have discussed in this paper. 
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In particular, supernovae type Ib,c are the result of core collapse of relatively small stars, 
where the hydrogen and perhaps also the helium shells are missing. A typical radius for 
SN Ic is R ~ 10 4 km. Therefore, the axion decay length D a = c T7 is comparable to R. 
This in turn implies that exp(—R/D a ) is not very small, it can be of order 0.1 to 1, and 
thus a reasonable amount of axions decay outside the mantle. Clearly, this also implies 
that the bulk of axions decaying outside the star produces a fireball (with a Lorentz factor 
T ~ 3 T/2 m e ). We can conclude that, associated with a supernova type Ib,c a GRB 
could take place, having a typical energy ~ 10 51 erg. 

In case of SN type II, which are associated with large stars, e.g. red giants, having an 
extended hydrogen shell (R ~ (1 — 20) x 10 7 km), the axion decay takes place completely 
inside the mantle and thus cannot be observed as GRB. For example, for SN 1987A, 
related to the blue giant Sanduleak having R = 3 • 10 7 km, the fraction of axions decaying 
outside the envelope, exp(— R/D a ), is essentially zero. Nevertheless, it can be of interest 
to think that the energy of order 10 51 erg released after axion decay at distances of few 
hundred or thousand km can help to solve the painful problem of mantle ejection and 
therefore of supernova explosion (in the prompt mechanism, shock stalls at a distance of 
few hundred km).Q 



6. Summary and conclusions 

We have shown that it exists a range of values for the PGB parameters such that a strong 
PGB emission takes place during the NS-NS or NS-BH merger. Choosing appropriately 
the parameters' value, the total luminosity for the PGB emission can be as large as few x 
10 55 erg/s. Moreover, at variance with the mechanism based on vv — > e + e~ annihilation 
which has a small (about one percent) efficiency, the PGB emission mechanism suffers no 
energy deficit: all emitted PGBs, due to their decay into e + e~ or photons, convert into 
the ultrarelativistic plasma, the fireball, which is at the origin of the burst. This plasma 
can be produced at distances as large as 100 km < D a < 10000 km (which are difficult 
to achieve in the case of plasma produced via vv annihilation), is not contaminated by 
baryonic matter and thus a large Lorentz factor can be obtained. Our conclusion is 



supported by the numerical simulations given in the recent ref. ||15|| . 

An interesting possibility opened by our mechanism is the association of "weak" GRBs 
with the supernovae type Ib,c. In addition, the PGB emission could also help to solve 
the problem of supernovae type II explosion, thus providing an unified theoretical base 
for the GRB and SN phenomena. 

The needed parameter range for such a PGB indicate an axion-like particle having a 
mass m a of a few MeV, coupling to nucleons g a N ~ few x 10~ 6 and to electrons g ae ~ 
few x 10~ 9 . This range of coupling constants coincides to that of the invisible axion with 
the Peccei-Quinn symmetry breaking scale / ~ few x 10 5 GeV. However, for such scale / 
a "true" axion would have a mass of few eV whereas we our PGB has a mass of few MeV, 



6 The possibility that the two photon decay of the axion-like particles with mass 0.15 — 1 MeV could 
provide a mechanism whereby gravitationally collapsing massive stars may eject their outer mantles and 
envelopes in supernova explosions was hrst considered in ref. |31|. It should be noted, however, that the 
range for the coupling g ai explored in ref. pl| | actually corresponds to the vertical band of our Fig. 3 and 
it is strongly excluded by the experimental limits |29|] . We thank S. Blinnikov for driving our attention 
to the ref. 
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i.e. about a million times larger. Present experimental data and astrophysical limits 
cannot exclude this particle, however the relevant parameters' window is not far from the 
present experimental limits and it can be experimentally tested in the close future. 
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